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Abstract 

We examine the optical properties of metals containing a periodic arrangement of 
nonoverlapping spherical mesopores, empty or filled with a dielectric material. We 
show that a slab of such a porous metal transmits light over regions of frequency 
determined by the dielectric constant of the cavities and the fractional volume 
occupied by them, with an efficiency which is many orders of magnitude higher than 
predicted by standard aperture theory. Also, the system absorbs light efficiently over 
the said regions of frequency unlike the homogeneous metal. 
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Recently, optical transmission experiments [1] have shown that periodic ar- 
rays of cylindrical holes in metallic films display a transmission efficiency, at 
wavelengths much larger than the hole diameter, which is orders of magnitude 
higher than predicted by standard aperture theory [2] . Such intriguing optical 
properties may have important technological applications in photolithogra- 
phy, in near-field microscopy, in flat-panel displays, and in novel active filter 
devices [1,3]. Similar enhanced optical transmission occurs in metallic grat- 
ings with very narrow slits [4], and it has been recognized that the physical 
mechanism behind this phenomenon relies on resonant transmission through 
coupled surface-plasmon modes. 

It has long been known that when metals are bombarded with energetic par- 
ticles such as neutrons or ions over a sufficiently long time, regular arrange- 
ments of vacancy clusters may form, such as three-dimensional (3D) lattices of 
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spherical voids or bubbles. The lattice of voids or bubbles is isomorphic to the 
microscopic lattice of the metallic matrix, with a lattice constant which is typ- 
ically 100 to 1000 times larger than that of the matrix lattice (see, e.g., Ref. [5] 
and references therein). Nowadays, considerable advances are being made in 
the template-assisted assembly of macroporous and mesoporous metals: a 3D 
template is assembled from a self-organizing material and impregnated with 
the desired metal; then, the template is removed, resulting in an array of 
pores that reflects the structure of the template [6]. In this paper we predict 
some extraordinary optical transmission effects in such mesoporous metals. 
We show that these properties can be understood in a systematic manner on 
the basis of exact calculations, which at the end provide a transparent model 
of the underlying physics. 

We consider, to begin with, a single spherical cavity, of radius S, in a metal 
characterized by a Drude-like relative dielectric function 

6=1- UjI/uj' , (1) 



where Up is the bulk plasma frequency, and we have neglected damping for 
now. The electromagnetic (EM) field at frequency uu is described by its electric- 
field component E(r; t) = Re[E(r) exp{—iijt)], where E(r) is written as follows 
(see, e.g., Ref. [7]). Inside the cavity (r < S) 
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where Kg = uj/c, c being the velocity of light in vacuum; ji is a spherical Bessel 
function; and X;^(f) is a vector spherical harmonic. We need not write down 
the associated magnetic-field component of the EM wave. Outside the cavity 
{r>S) 
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where hf' is a spherical Hankel function corresponding to an outgoing wave. 
The first two terms in the above equation describe an incident wave and 
the last two terms a scattered wave. The wavenumber k for a medium with 
negative dielectric function, as is the case to be considered here {cu < cUp), is a 
purely imaginary number: k — iq — icuy/^/c. 



Because of the spherical symmetry of the scatterer we obtain 
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where 
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with a corresponding expression for T^^ {uj) [7]. In the case of a single cavity in 
a homogeneous medium of negative dielectric function there can be no incident 
wave: = = and, therefore, non-trivial states of the EM field, of given 
I, will exist if 

Ti^{u) = oo or Ti"{u) = oo . (6) 



It can be shown that the second of the above equations does not obtain a 

real frequency solution in the region from to Up. The first equation, for any 
/, has a solution in this frequency region; we denote it by Coi. For a small 
cavity, using the appropriate asymptotic expansions of the spherical Bessel 
and Hankel functions for small arguments [8] in Eq. (5), we obtain 

(I;,~u;,^(Z + l)/(2Z + l), Z = l,2,3,.... (7) 



Taking for instance Sup/c — 0.4, which for cUp — 10 eV corresponds to a 

radius 5" = 7.5 nm, Eq. (7) approximates ui with an accuracy better than 
1%. The eigenmodcs of the EM field obtained at a)/ define the 2^-pole plasma 
oscillations at the surface of the cavity. 

Let us now ask whether it is possible for a system of small spherical cavities in 
a metal to be transparent. For a crude, qualitative description of the optical 
properties of the composite system, we can use its effective dielectric constant, 
e, as given by the Maxwell Garnett (MG) effective-medium theory [9] 

^ = /^^, (8) 
e + 2e + 2e ' ^ ^ 



where / is the fractional volume occupied by the cavities (e^ = 1) and e is 
the relative dielectric function of the metal given by Eq. (1). Eq. (8) gives a 

positive e over the range of frequencies from dJmin = c<jpy^2(l — f )/3 to cj^ax = 

cUp^J (2 + f)/3. In other words, in the above frequency region, there is a band 
of propagating states of the EM field in the porous metal, the width of which, 

<^max - <^min, inCrCaSCS with /. As / ^ 0, Umax ^ <^mm ^ ^p\/V^ which. 
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according to Eq. (7), is the resonance frequency, uji, associated with the dipole 
plasma oscillations at the surface of a single small cavity. Clearly the above 
band (for / > 0) arises from the interaction between dipole plasma modes 
localized on neighboring cavities and can therefore be understood in the spirit 
of the tight-binding approximation [10]. Its bandwidth will be larger the larger 
the spatial overlap between the wavefields associated with neighboring spheres 
which in turn increases with /. 

The entire band can be shifted to lower frequencies by filling the pores with 
a dielectric material. One can show that the resonance frequencies of a single 
small sphere of dielectric constant in a metal described by Eq. (1) are given 



Therefore, with ~ 10, Ui of the isolated sphere and the corresponding band 
of the porous metal are found in the optical region. 

We note that the MG theory is based on the electric-dipole approximation 
and, therefore, cannot describe 2'-pole states with I > 1. There is of course 
an infinite number of bands corresponding to the 2'-pole plasma modes of the 
individual spheres. These modes, according to Eq. (9), extend approximately 

from UpJl/{es + 1) to ujpJ2/{es + 2). It is therefore clear that there are no 



frequency bands below a cutoff frequency at about cj^y 1/(6^ + 1). This is what 
one expects of a metal (a system with DC conductivity at any temperature in- 
cluding zero), and the system under consideration possesses DC conductivity 
because the metallic component of it forms a continuous network [11]. It is un- 
derstood, however, that at low frequencies {lot < 1) Eq. (1) must be replaced 
by Eq. (10) below, which takes into account the damping of conduction-band 
electrons. We note that, in a system consisting of nonoverlapping metallic 
spheres in a dielectric medium, there are propagating modes of the EM field 
at low frequencies with a free-photon- like dispersion [12], reflecting the fact 
that the system is an insulator at a; = 0. 

An accurate analysis of the optical properties of porous metals can be effi- 
ciently carried out using the method we developed for the calculation of the 
optical properties of photonic crystals consisting of nonoverlapping spheres in 
a homogeneous host medium. The details of the method and a computer pro- 
gram for its implementation can be found elsewhere [13,14]. Here we need only 
say that in the present case, of a metallic host medium, where only evanescent 
EM waves exist below the plasma frequency, the f2-matrices (as defined in 
Ref. [14]) must be calculated by a direct summation over the space lattice; a 
few terms in the relevant lattice sums give adequate convergence. 
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Wc consider a system of nonoverlapping, identical silicon spheres (e^ = 11.9), 
centered at the sites of a fee lattice of lattice constant a, in a metallic host 
medium whose dielectric function is given by Eq. (1). The radius, S, of the 
spheres equals one fifth of the first- neighbor distance, Oq, which corresponds to 
a fractional volume occupied by the spheres / = 4.74%. We view the crystal 
as a stack of layers (planes of spheres) parallel to the fee (001) surface. For 
given k||, the reduced component of the wavevector parallel to the fee (001) 
surface, we calculate, as functions of cu, the frequency lines — /c^(cj;k||) 
corresponding to generalized Bloch wave solutions of the EM field. kz{uj; ky) 
is the z component [normal to the (001) plane] of the wave vector of a gen- 
eralized Bloch wave with the given uj and ky; and there are many such waves 
corresponding to different values of kz- The regions of uj over which at least 
one kz is real define corresponding frequency bands and regions over which all 
kz are complex define frequency gaps, for the given k||. 

In Fig. la we show the frequency bands for ky = [dispersion curves along the 
normal to the (001) plane] of the above crystal, which arise from the 2'-pole 
states of the individual spheres for I = 1,2, 3. We see that these bands develop 
about the corresponding resonance frequencies given by Eq. (9), and that their 
width decreases with I. Apparently the spatial extent of the wavefield around 
a sphere decreases with increasing thus leading to a weaker hybridization 
between the higher 2'-pole states of the spheres and, consequently, to a smaller 
bandwidth. The bands shown by the black solid lines in Fig. la are doubly 
degenerate and couple with light incident normally on a slab of the crystal par- 
allel to the (001) surface. The black broken fine shows the dipole band (doubly 
degenerate) as obtained in the MG approximation. The bands shown by the 
gray solid lines are nondegenerate and do not couple with normally incident 
light. The existence of these optically inactive modes at high-symmetry points 
and along symmetry lines has been noted by a number of authors [13,15] . Fig. 2 
shows more clearly all frequency bands corresponding to I — 2 and 1 — 3. One 
can see that the total number of bands associated with a given / equals 21 + 1 
as expected from the {21 + l)-degeneracy of the corresponding state of the 
isolated sphere. The degeneracy of the bands and the symmetry properties of 
the corresponding eigenmodes is what one expects from a group-theoretical 
analysis [16]. Next to the frequency band structure, in Fig. lb, we show by 
the solid line the transmission coefficient of light incident normally on a slab 
of the crystal consisting of 16 planes of spheres parallel to the (001) surface, 
assuming that the medium on either side of the slab is air. As expected, the 
transmittance practically vanishes for frequencies within the frequency gaps of 
the infinite crystal and exhibits the well-known Fabry-Perot-type resonances 
over the regions of the optically active frequency bands [12] . The broken line 
shows the transmittance as evaluated in the MG approximation. 

The absorption of light by the material, which we have neglected so far, is 
of course a quantity of great importance. We account for it by inserting, as 



5 



usual, a damping term in the relative dielectric function given by Eq. (1). We 
write 
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where r is the relaxation time of the conduction-band electrons, with a typical 
value {u!pT)~^ = 0.001. Absorption is mainly associated with the resonantly 
oscillating multipoles of the spheres and, therefore, is appreciable only in the 
region of the frequency bands. This is demonstrated in Fig. 3a, for light in- 
cident normally on a slab of the crystal consisting of 16 planes of spheres 
parallel to the (001) surface. Wc sec that absorption is considerable; the cor- 
responding absorbance of a homogeneous metallic slab of the same thickness 
is of the order of 10"^. In Fig. 3b we show the transmission coefficient of light 
for the same system; transmission occurs in the frequency region of the dipole 
band. It is somewhat reduced, due to absorption, but it is extraordinarily high: 
the corresponding transmittance of a homogeneous metallic slab of the same 
thickness is of the order of 10~^, while a single hole of diameter o? in a metal 
film transmits light of wavelengths X>> d with an efficiency proportional to 
{d/Xy [2] and in our case we have d/X^ 0.02. 

In addition to the above, the results shown in Figs. 1 and 3 indicate that 
though the MG approximation gives a reasonable estimate of the transmit- 
tance and absorbance of the system in the frequency region of the dipole 
band, it differs significantly from the exact picture. And of course, it fails to 
reproduce all features associated with the higher 2'-pole resonances (see also 
Ref. [12]). Finally, it appears that the position of the frequency bands (and 
consequently the range of frequencies over which a high transmission and ab- 
sorption occurs) depends mostly on the characteristics of a single sphere in 
the metal host and on the fractional volume occupied by the spheres, and less 
so on the specific arrangement of the spheres. We should also emphasize that 
though our results have been obtained for a metallic component described by 
Eqs. (1) and (10), what really matters is that the dielectric function of the 
host medium is such that solutions of Eqs. (6) exist at appropriate frequencies; 
and as long as this remains true, similar results will be obtained. 
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Fig. 1. (a) The photonic band structure normal to the (001) surface of a fee crystal of 
silicon spheres in a nonabsorbing Drude metal (aQUJp/c = 1, Su>p/c = 0.2). The solid 
lines (black lines: doubly degenerate bands, gray lines: nondegenerate bands) are 
exact results; the broken line is the MG result, (b) The corresponding transmittance 
curve for light incident normally on a slab of the above crystal consisting of 16 lattice 
planes parallel to the (001) surface. The solid (broken) lines are exact (MG) results. 




Fig. 2. The frequency bands associated with the Z = 2 (a) and Z = 3 (b) states of 
an isolated sphere, for the system described in Fig. 1. 
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Fig. 3. Absorbance (a) and transmittance (b) of light incident normally on a slab 
of 16 lattice planes parallel to the (001) surface of a fee crystal consisting of silicon 
spheres in a Drude metal [a^Up/c = 1, Sujp/c = 0.2, [ujpT)~^ = 0.001]. The solid 
(broken) lines are exact (MG) results. 
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